Moment problems and boundaries of number triangles 



oo 

o 



^1- 



Alexander Gnedin* and Jim Pitman''' 
February 22, 2008 



Abstract 



The boundary problem for graphs like Pascal's but with general multiplicities of edges is related 
to a 'backward' problem of moments of the Hausdorff type. 
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1 The extreme boundary 

Let T n := {(n, 0), (n, 1), . . . , (n,n)} and T := Li^L T n . We endow T with the structure of a directed 
Oh ' graph in which every node (n, k) has two outgoing edges (n, k) — ► (n + 1, k) and (n,k) — > (n + 1, fc + 1) 

with multiplicities £ n k and r n k (respectively), where {£ n k] {n,k) £ T} and {r n jt; (n, k) £ T} are given 
triangular arrays with (strictly) positive entries. A classical example is the Pascal graph with unit 
. multiplicities £ n k = r n k = 1. 

i *~ l i ' Let V be the set of nonnegative solutions V = {V n k'i (n, k) £ T} to the backward recursion 

7— I ■ Kfc = (-nkVn+i.k + r nk V n+ i tk , {n,k)eT (1) 

> ; 

with normalisation Voo = L The set V is convex and compact in the product topology of functions on 
T . By some general theory in Dynkin (1978) the extreme boundary extV, comprised of indecomposable 
elements of V, is a Borel set. Moreover, V is a Choquet simplex, meaning that each V S V has a unique 
representation as convex combination 

O I V= I U^{dU) (2) 

QQ ' J CXtV 

with some probability measure /i supported by extV. The boundary problem for the graph T is to find 
some explicit description of the set of extremes, meaning, if possible, a simple parametrisation of extV 
• ^ • along with the kernel that is implicit in Q . 

The recursion |T]) for the Pascal graph appeared in the work of Hausdorff on summation methods (1921, 
p. 78) and the 'small' problem of moments on [0, 1], In this case the bivariate array V is completely 
determined by V,.o according to the rule V,+k,k = V fe (V #i o), where • stands for the variable n, and V fc is 
the fcth iterate of the difference operator V(E/.) := U, — C/.+i. The condition V > means that V,.q is 
completely monotone, hence by Hausdorff's theorem V,fi is a sequence of moments 



V n0 = / x n fi(dx) 

J[0,1] 

of some probability measure [i. That is to say, the set of extremes extV can be identified with the unit 
interval, and the extremes have the form V n k(x) — x n ~ (1 — x) for x £ [0, 1]; in particular, V n o — x n . 
For general multiplicities the recursion JTJ) is equivalent to 

V.,k = Vfc(---(Vi(V , . |0 ))---), 
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where Vfe(J7.) = (U, — £,,kU,+i)/r,,k is a generalised difference operator. By analogy with Hausdorff's 
criterion, the question about positivity of the generalised iterated differences of V #i o may be regarded as 
a 'backward' problem of moments. A direct problem of moments of the Hausdorff type appears when 
we determine the V n o's for extreme solutions as functions on the boundary, and consider the integral 
representation of the generic V n0 in the form @ . 

A bivariate array V € V could be also computed by suitable differencing the diagonal sequence (V nn ), 
but this leads to the same type of the moment problem by virtue of the transposition of T which exchanges 
the multiplicities £ n k's with r„ iT1 _fe's. 

A special feature of T, as compared with more complicated graphs like Young's lattice (see Kerov 
(2003), Borodin and Olshanski (2000)), is a natural total order on the extreme boundary. In this note 
we extend the argument of Gnedin and Pitman (2006) to show that the total order allows extV to be 
embedded into [0,1]. We shall also survey the connection of the boundary problem with asymptotic 
properties of some classical arrays of combinatorial numbers. 

2 Markov chain approach 

The weight of a path in T joining the root (0, 0) and some other node (n, k) is defined as the product of 
multiplicities of edges along the path (for instance the weight of (0,0) —> (1,0) — > (2, 1) is A)o r io)- The 
dimension D n f, of (n, k) £ T is defined to be the sum of weights of all paths from (0, 0) to (n, k). The 
dimensions are computable from the forward recursion 

D n k = r n -i,k-i Ai-i,fc-i + t n -i,kD n -i,k, ( n : k) £ T, (3) 

(where the first term in the right-hand side is absent for k — and the second term is absent for k = n), 
with the initial condition Dqq = 1. The number triangle associated with T is the array {D n k; (n, k) E T}. 

Each V 6 V determines the law Py of a inhomogeneous Markov chain K, whose backward transition 
probabilities for < fc < n, n > are 

P v {K n -i = j | K n =k) = + r n _i,A*_i), (4) 

and whose distribution at time n is ¥y{K n = k) = Dn~kVnk- It is important that the probabilities ^ are 
determined solely by the multiplicities of edges and do not depend on V. Hence V is in essence a class of 
distributions for Markov chains on T with given backward transition probabilities. 

For each fixed integer v and < x < v let V V3t be the function on T which satisfies the recursion {T]) 
for n < v, satisfies V"j* = 6^, and equals on \J n>v T n . Such V V7< determines the probability law of a 
finite Markov chain (Kq, . . . , K v ) conditioned on K v = x. 

We define the sequential boundary d-\V to be the set of elements of V representable as limits V — 
lim^oo V"'*^") taken along infinite paths \k(v)\ v = 0, 1, . . .} in T. The sequential boundary <9jV may 
be smaller than the set of all accumulation points for {V v>c \ [y, x) G T} (the Martin boundary), but it 
is large enough to cover extV, as is seen from the following lemma, which is a variation on the theme of 
sufficiency (see Diaconis and Freedman (1984)). 

Lemma 1 . IfVG extV then the random functions 

v»> K » : = ]T i {Jt ^ } r (5) 

3<=0 

satisfy V V ' K " — > V, as v — > oo ; Py -almost surely. 

Proof. Let T v be the sigma-algebra generated by K v , K v +i, . . ., and = V\T V . Let Py correspond to 
some extreme V. Choose any (n, k) and consider random variables 

Kf" = ^v{K n = (n, k) | K v )/D nk = ¥ v (K n = (n, k) \ F u )/D nk , v>n 
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where the first equality follows from the definition ((5]), and the second equality is a consequence of 
the Markov property. Applying Doob's reversed martingale convergence theorem to the conditional 
expectations given f i D J2 3 . . . we obtain 

Kf" -» IV {K n = (n, k) I Foo)/D nk P y -a.s. 

The assumption V G extV implies that is trivial, hence 

Py(X„ = (n,k) |^oo) = Py(X„ - (n,fe)) - 



Thus extV C 9|V (in general the inclusion is strict). To state this conclusion in analytical terms, 
define the weight of a path in T connecting two nodes (n, k) and (u, x) as the product of multiplicities 
along the path, and define the extended dimension as the sum of weights over all such paths (so that 
Dqq = D v>c ). We then have a fundamental relation 

= (6) 

which connects the boundary problem with asymptotic properties of T. Specifically, the convergence of 
yv,K(u) amoun t s t the convergence of these ratios for all (n, k) G T along the path (in fact, it is enough 
to focus on K,o)- 



3 Order 

A special feature of T which yields the order is that the only possible increments of the variable k along 
any path are and 1. The next lemma appeared in Gnedin and Pitman (2006) with a different proof. 

Lemma 2. For v > n fixed, V^q is nonincreasing in x. 

Proof. Choose < x < x 1 < v and consider two Markov chains K m , K' m which run in reverse time 
n = v, v — 1, ... ,0 according to ^ and start with K v — x,K' v = xf . Suppose the chain K' m jumps 
independently of K, as long as they are in distinct states, and suppose that K' m is coupled with K m at 
some random time < r < v when the states become the same. In the reverse time, only transitions 
k — > k, k — ► k — 1 for k > and 0^0 are possible, hence we always have K' n > K n . Therefore the event 
K' n = occurs exactly when K n — and r > n, which implies 

F{K n =0\K v = x)> ¥(K n = 0\K U = x'). 



A minor modification of the above argument shows that if K v under Py is strictly stochastically 
smaller than K v under some other Py, then the same relation holds true for every n < v. 

We focus now on Vi . Suppose V E d^V is induced, via ((6]), by some infinite path {x{v); v = 0, 1, . . .}, 
and V' G 9fV is induced by some other path {x 1 {v); v = 0, 1, . . .}. If x(y) = x 1 (y) for infinitely many 
v then, of course, V — V. If x(y) < x'(y) for infinitely many v and x(y) > x 1 (y) for infinitely many 
v then by Lemma [2] we have V, t o = V.' and V = V . Thus V ^ V can only occur if the same strict 
inequality holds for all sufficiently large v. To be definite, let x(v) < xf{v) for all large enough v, but 
then V 7^ V implies that K n under Py is strictly stochastically smaller than K n under Py for all n > 0, 
in particular this holds for n = 1 which means that Viq > V{ Q . We see that for V, V G 9j V, the inequality 
Via > V( holds if and only if K n under Py is strictly stochastically smaller than K n under Py/ for all 
n > 0. This defines a strict order <a on 9|V. 

Lemma 3. The sequential boundary d^V is compact. 

Proof. Suppose V J G 9|V (j = 1,2,...) is a sequence converging to some V G V. We know that V 
is a metrisable compactum with some distance function dist. Passing to a subsequence we can restrict 
consideration to the case of increasing or decreasing sequence, so to be definite assume that V^ +1 <\ V-? 
for j — 1,2,... Choosing some path {x^v); v = 0, 1, . . .} which induces V J , the ordering implies that 
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x-- 7 (u) — > oo as v — > 00 and x^{y) < x^ +1 (v) for all large enough ^. As ^ varies, define inductively in j a 
function x(f) which coincides for some v with >f J '(z/). Specifically, x(y) = x^iy) until x^ +1 {v) < x^(y) 
starts to hold along with d\st{V u,><i , V) < 1/j and dist(y ,>fj+1 , V) < l/j, then let x(v) decrement by 1 
until it becomes equal to x J+ (y). This defines an infinite path in T, for which one can use monotonicity 
to show that V V '"M -> V. 

Recalling that ^oo^io + r oo^ii = 1 we obtain: 

Theorem 4. The function V i— ► £00^10 is an ordered homeomorphism of the sequential boundary dfV 
with order < into [0, 1] wrai/i order >. 

Two extreme cases £qqVq\ = and £00 7)1 = 1 correspond to trivial Markov chains K, = and K, = •, 
respectively. 

4 Discrete or continuous? 

In the situation covered by the following lemma, setting x{v) = m (for large v) for m = 1, 2, ... is the 
only way to induce nontrivial limits. Then extV is discrete and coincides with the sequential boundary. 

Lemma 5. (Gnedin and Pitman (2006)) Suppose for m = 0, 1, . . . there are solutions V(m) € V such 
that V nm (m) D nm — ► 1 as n — > 00, t/ien eac/i V(m) is extreme and satisfies K n — > to Py( m )-a.s.. J/ also 
Vio(»7i) — * tts to — > 00 i/ien V(m) converges to the trivial solution V(oo) with K, = • Py( oc )-a.s., and 
in this case extV = <9 T V = {V(0), 7(1), . . . , V(oo)}. 

In some cases the limits can be obtained by setting x[y) ~ s c(v) with suitable scaling c{v) — > 00 and 
s > 0. Under conditions in the next lemma, extV coincides with c^V and is homeomorphic to [0, 1]. The 
scaling determines the order of growth of K, under Py's. 

Lemma 6. (Gnedin and Pitman (2006)) Suppose there is a sequence of positive constants {c(v)]v = 
0, 1, . . .} with c(y) — > 00, and for each s £ [0, 00] there is a solution V(s) S V which satisfies K v jc{y) — > s 
Py( s )-a.s.. Suppose the mapping s 1— > V(s) js a continuous injection from [0, 00] to V im£/i and 00 
corresponding to the trivial solutions. Then a path {x(i>); v = 0,1,...} induces a limit if and only if 
x(y)jc(y) — > s for some s G [0, 00], in which case the limit is V(s). Moreover, extV = dV\ = {V(s), s £ 
[0,oo]}. 

Minor variations of the above two situations are obtained by transposing multiplicities £ 
Still, this does not exhaust all possibilities. See Kerov (2003) (Section 1.3, Theorem 2) for examples of 
boundaries with both discrete and continuous components. 

5 Number triangles 

The Pascal triangle. For the Pascal graph the dimensions are D n k = (^) and = (^Zt)- The 

ratios V^j. = (j(uy^k)/(>c[v)) conver g e iff *c{y)jv — > x € [0,1], in which case the limit is V n o(x) = x n . 
This identification of extremes is equivalent to de Finetti's theorem (see Aldous (2003)), since V £ V 
determines the law of some infinite sequence of exchangeable Bernoulli trials. A closely related type of 
moment problem with a monotonicity constraint have been discussed recently in Gnedin and Pitman 
(2007). 

The o-Pascal triangle. This graph has multiplicities £ n k — l,r n k = q n ~ k , (n,k) G T, and may be 
seen as a parametric deformation of the Pascal graph. The extreme boundary was found in Kerov (2003) 
by an algebraic method and justified by Olshanski (2001) by the analysis of ©■ The dimensions are 
expressible through g-binomial coefficients as 

Suppose first that < q < 1. Lemma [5] is applicable, and all nontrivial extremes are given by 

g("-fc)("-fc)(l-g)...(l-g") 

V nk (m) = — — — — — — ^ l {0 <fe<m}, m = l,2,... 
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In particular, V n o(m) = q mn , (m = 0, 1, . . . , oo). 

The function V i— > ^oo^io identifies the extreme boundary with {q m , m = 0, 1, . . . , oo}. The decom- 
position © into extremes corresponds to a version of Hausdorff 's moment problem on [0, 1] with kernel 
x n , but subject to the constraint that the measure is to be supported by a geometric progression. That 
is to say, a sequence V, o with Voo = 1 is representable as a mixture 

t t \ ^ m mn 

m£{0,l,...,oo} 

with some probability distribution {p m ; Tn — 0, 1, . . . , oo} if and only if = Vfe(- • • (Vi(V, ! o) • ■ • ) > 
for all k > 0, where V fe (f7.) = (Z7. - J7. + i)/g'^ fc . 

In the case q > 1 the extreme boundary is {1 — q~ m , m = 0, 1, . . . , oo} (this case is reducible to q < 1 
by transposition of T and replacing q by g -1 ). The only accumulation point of extV for q < 1 is and 
for q > 1 is 1. A phase transition occurs at q = 1, when the extreme boundary is continuous. 
Stirling triangles. Let r n f. = 1 and £„fe = ( n + 1) — + 1) f° r — 00 < a < 1. For a = — oo take 

n + 1 1 

The notation stands for the generalised Stirling numbers 



Axfc = fc + 1- The dimension is D n k = 
defined as connection coefficients in 



k + 1 



fc=i 



l (t/a) nh 



(where f denotes the rising factorial), with the convention that these are the Stirling numbers of the 
second kind for a — — oo. For a — these are the signless Stirling numbers of the first kind. 
For — oo < a < the extreme boundary is discrete, with 

V n0 (m) = - — 1 for - oo < a < 0, V n o(m) = -\ for a = -oo. 

[m\a\ + l) n m" 

These kernels under ly a moment problem for measures on the set {0,1, ...,oo}. 

A phase transition occurs at a = 0. Lemma [5] applies with x{v) ~ s logn, the extreme boundary is 
continuous and the kernel is 

Ko(s) = , tt s e [0, oo]. 

(s + 1)„ T 

This case is closely related to random permutations, records and Ewens' sampling formula (see Arratia 
ct al (2003)). 

In the case < a < 1 we should take x(v) ~ sn a to generate the boundary, see Gnedin and Pitman 
(2006) for formulas for V n o(s) (to adjust the notation in Gnedin and Pitman (2006) to the present setting, 
one should replace (n, k) by (n + l,k+ 1)). This family of solutions is related to Poisson-Kingman 
partitions, see Gnedin and Pitman (2006) and references therein. 

Several results and (still open) conjectures about boundaries of more general Stirling graphs, with 
multiplicities of the form £ n k = b n + a./., r n y. = 1, are given in Kerov (2003). 

The Eulerian triangle. For multiplicities £„fe = k + 1, r n j~ = n — k + 1 the dimension is the Eulerian 
number i"^ 1 ) (that counts permutations with a given number of descents). The boundary problem was 
solved in Gnedin and Olshanski (2006). The extreme solutions are given by 

^ n — k / 

V n k(m) = - — TTTu II + - 
[n + 1)! - 1 -■- V m 

v i— — k v 

with m£ZU {oo}. Note that the range of £ooVio(m) = (m + l)/(2m) is symmetric about 1/2, with 1/2 
being the only accumulation point. The symmetry of the boundary stems in this case from the invariance 
of multiplicities under transposition. 
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